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The production and propagation of antikaons - described by dynamical spectral 
^ ' functions Afi{X, P, M'^) as evaluated from a coupled channel G-matrix approach - 

C ■ is studied for nucleus-nucleus collisions at SIS energies in comparison to the con- 

ventional quasi-particle limit and the available experimental data using off-shell 
transport theory. We find that the K~ spectra for ^'^C+^'^C and ^^Ni+^^ Ni at 1.8 
A-GeV remain underestimated in the G-matrix approach as in the on-shell quasi- 
particle approximation whereas the preliminary spectra for Au + Au at 1.5 A-GeV 
are well described in both limits. This also holds approximately for the K~ rapidity 
distributions in semi-central collisions of Ni + Ni at 1.93 A-GeV. However, in all 
limits considered there is no convincing description of all spectra simultaneously. 
Our off-shell transport calculations, furthermore, demonstrate that the strongest 
in-medium effects should be found for low antikaon momenta in the center-of-mass 
frame, since the deceleration of the antikaons in the attractive Coulomb and nu- 
clear potentials and the propagation to the on-shell mass induces a net shift and 
squeezing of the K~ spectra to the low momentum regime. 
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1 Introduction 



An open problem of todays strong interaction physics is the dynamical generation of 
hadron masses from the quark and gluon fields, that are the elementary fields of Quantum- 
Chromo-Dynamics (QCD). Whereas QCD lattice calculations nowadays give quite reliable 
results for the masses of the 0~ and 1~ meson octet states - except for the light pions 
- as well for the baryon octet and decuplet states, the masses and properties of these 
states are not well known at finite quark chemical potential /ig or finite baryon density 
Pb- This is essentially due to the problem of calculating the fermion determinant at finite 
Hq. Though some recent progress has been made for QCD lattice calculations at finite 
baryon density [1, 2], the properties of hadrons in a baryonic environment cannot be 
addressed directly so far by such ab-initio calculations. One thus presently has to rely 
on effective hadronic Lagrangians and to compute the hadron properties at finite density 
Pb and/or temperature T by 'dressing' the vacuum states in the medium e.g. within 
G-matrix theory. An alternative way is to employ effective chiral Lagrangians and to 
extract the leading order effects for the hadron modifications in the medium. 

As demonstrated in the pioneering work of Kaplan and Nelson [3, 4] within a chiral 
approach kaons and antikaons couple attractively to the scalar nucleon density with a 
strength proportional to the KN — E constant, 

Sktv = +m°) < N\uu + ss\N >, (1) 

which may vary from 270-450 MeV (cf. the discussion in Rcf. [5]). QCD lattice calcula- 
tions here provide further information [6]. In Eq. (1) and m° denote the bare masses 
for the light u- and strange s-quark while the expression in the brackets is the expectation 
value of the scalar fight and strange quark condensate for the nucleon [7] . Furthermore, 
a vector coupling to the quark 4-current - for vanishing spatial components - leads to a 
repulsive potential term for the kaons; this (Weinberg- Tomozawa) term is attractive for 
the antikaons. 

In relativistic mean-field models the dispersion relation for kaons and antikaons in the 
nuclear medium - for low momenta - can be written as [8, 9, 10] 



U^K^pN.V) = ±075- + 

8 fi 



2 , 2/1 ^KN , ( ^PN 



1/2 

(2) 



In Eq. (2) rriK denotes the bare kaon mass, ~ 93 MeV is the pion decay constant, 
while ps and pN stand for the scalar and vector nucleon densities, respectively. Note that, 
when extrapolating Eq. (2) to 3po and above, the antikaon mass becomes very light. 
Thus antikaon condensates might occur at high baryon density which, furthermore, are 
of interest in the astrophysical context [11, 12, 13, 14]. 

Studies on production from nucleus-nucleus coUisions at SIS energies of 1- 2 A- 
GeV have shown that in- medium properties of kaons are seen in the collective fiow pattern 
of mesons, both in-plane and out-of-plane, as well as in the abundancy and spectra 
of antikaons [7, 15, 16, 17, 18, 19]. Thus in- medium modifications of the mesons have 
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become a topic of substantial interest in the last decade triggered in part by the early 
suggestion of Brown and Rho [20], that the modifications of hadron masses should scale 
with the scalar quark condensate < qq > at finite baryon density. 

The actual kaon and antikaon self energies (or potentials) are quite a matter of debate 
and depend on the many-body resummation scheme adopted. Especially for the antikaons 
the momentum-dependence of their self energies is widely unknown since most Lagrangian 
models restrict to s-wave interactions or only include additional p- waves [21, 22, 23]. Here 
only limited information is available from a dispersion analysis in Ref. [24] . 

There have been attempts to extract the antikaon-nucleus potential from the analysis 
of kaonic-atom data and some solutions favor very strongly attractive potentials of the 
order of -200 MeV at normal nuclear matter density po [25, 26]. However, more recent self- 
consistent calculations based on a chiral Lagrangian [27, 28] or meson-exchange potentials 
[29] only predict moderate attractive depths of -50 to -80 MeV at density po- In addition, 
studies of kaonic atoms using the chiral KN amplitudes of Ref. [30] show that it is indeed 
possible to find a reasonable reproduction of the data with a relatively shallow antikaon- 
nucleus potential [31, 32], albeit adding an additional moderate phenomenological part 
[33]. This has been corroborated by a calculation in Ref. [34], where a good fit to both 
scattering K~p data and kaonic-atom data only required to modify slightly the parameters 
of the chiral meson-baryon interaction model of Ref. [35]. Since it is clear that kaonic 
atom data do not suitably constrain the antikaon-nucleus potential at normal nuclear 
matter density, a recent work [36] explored the possibility of distinguishing between deep 
or shallow potentials [25, 28, 33] by means of nuclear scattering of low energy K~ produced 
from $ decay at Da$ne. 

In fact, the antikaon-nucleon amplitude in the isospin channel / = is dominated 
by the A (1405) resonant structure, which in free space appears only 27 MeV below the 
KN threshold. This resonance is generated dynamically from a coupled channel T-matrix 
scattering equation using a suitable meson-baryon potential. The coupling between the 
KN and ttY {Y — A, E) channels is essential to get the right dynamical behavior in 
free space. Correspondingly, the in-medium properties of the A(1405), such as its pole 
position and its width, which in turn infiuence strongly the behavior of the antikaon- 
nucleus optical potential, are very sensitive to the many-body treatment of the medium 
effects. Previous works have shown, for instance, that a self-consistent treatment of the 
K self energy has a strong influence on the scattering amplitudes [27, 28, 29, 37] and, 
consequently, on the in-medium properties of the antikaon. Moreover, the incorporation 
of the pion with its medium modified properties also proved to be an important aspect 
[28, 38], although most works until now have ignored it. As pointed out in Ref. [39] also 
the properties of the E*(1385) might be essential for in-medium transitions rates at finite 
relative momentum since the E*(1385) is the strange analogue to the A(1232). 

Furthermore, a sizeable progress in the understanding of hadron off-shell propagation 
in phase-space - even for complex configurations - has been obtained in the last years. In 
Refs. [40, 41] a semiclassical transport approach has been developed on the basis of the 
Kadanoff-Baym equations that includes the propagation of hadrons with dynamical spec- 
tral functions. This approach has been examined for nucleus-nucleus collisions at GANIL 
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[40], SIS and AGS energies [41] as well as for equilibration phenomena in related infinite 
nuclear matter problems [42]. Formally, the input for the off-shell transport approach 
with respect to the complex hadron self energies and the various channel transition prob- 
abilities should be provided by coupled-channel G-matrix calculations. Though being far 
from complete, the present work makes a major step in this direction by incorporating 
the off-shell information from KN ^ Yir transitions. 

The paper is organized as follows: In Section 2 we will briefly review the general- 
ized transport equations on the basis of the Kadanoff-Baym equations [43] and present 
the dynamical equations of motion for 'test-particle' propagation in 8-dimensional phase 
space in comparison to the traditional on-shell equations of motion. The input for this 
approach relates to self-consistent self energies for the hadrons which here are taken from 
the coupled-channel G- matrix approach following Ref. [38] (Section 3). The actual de- 
scription of in-medium kaon and antikaon production in NN and ttN collisions - in line 
with the K spectral function from the G-matrix approach - is described in Section 4. The 
explicit transport calculations for nucleus-nucleus collisions in the SIS energy regime are 
presented in Section 5 in comparison to the experimental data available. A summary and 
discussion of open problems concludes this study in Section 6. 

2 Description of the off-shell transport approach 

Whereas the equilibration of strongly interacting quantum systems has been studied on 
the basis of the Kadanoff-Baym equations [43] for infinite nuclear matter problems quite 
some time ago [44] the question how to propagate 'broad resonances' in an inhomogenuous 
medium out-of equilibrium (as encountered in relativistic nucleus- nucleus collisions), has 
been addressed and solved only in the last years [40, 41, 42, 45, 46, 47, 48, 49]. In 
the latter works a semi-classical off-shell transport approach has been derived from the 
Wigner transformed Kadanoff-Baym equations in the limit of first order gradients in 
phase space. The approach by Ivanov et al. in Refs. [46, 47, 48, 49] differs from those 
in Refs. [40, 41, 42, 45] in the treatment of the 'rearrangement term', which is kept in 
all orders in Refs. [46, 47, 48, 49], but only employs first order phase-space gradients in 
Refs. [40, 41, 42, 45] as suggested by Botermans and Malfliet already in 1990 [50]. The 
latter limit allows to perform a test-particle solution to the problem, which is adequate 
for present transport approaches, whereas the transport equations from Ivanov et al. 
[46, 47, 48, 49] might only be solved on a 8-dim. grid in phase space which is very 
unconvenient for actual applications. A recent overview on the different approaches has 
been given in Ref. [49]. 

We here follow the off-shell approach developed in Refs. [40, 41, 42]. For the actual 
details we refer the reader to the original articles or independently to Ref. [45] for a 
nonrelativistic formulation. We now concentrate on those results that arc important for 
our present study. First of all, the transport limit provides an algebraic result for the 
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hadron spectral function 



(P2 _ M2 - i?eS5,^*,)2 + r2,p/4 ' 

which holds for all approaches [40, 41, 42, 45, 46, 47, 48, 49]. In Eq. (3) Mq denotes the 
bare mass pole, Txp = — 2/mS5^p while S5^*p is the retarded self energy of the hadron 
which in general depends on space-time X and the four-momentum P in a hadronic 
environment. 



2.1 Test particle representation 



In order to obtain an approximate solution to the transport equation ((16) in Ref. [41]) 
we use a testparticle ansatz for the Green function Sxp, i-e. for the Wigner transform of 



7 S'< •- 



(4) 



with X = {x + y)/2, where ^{x),^''{y) denote the hadron field operators at space-time 
position x or y. More specifically, we rewrite iSxp in terms of the real and positive 
semidefinite quantity 



N 



S-xpNxp 



iS<p ~ Yl S^'\X - X,{t)) S^'\P - %)) S{Po- e,{t)) , (5) 



where the spectral function Axp [Eq. (3)] is separated from the number density function 
Nxp- In the most general case (where the self energies depend on four-momentum P, 
time t and the spatial coordinates X) the equations of motion for the testparticles read 
[41] 



dXj 
dt 



1 



1 — C(i) 2e 



2 P, + Vp, ReTTg + ^ '—-^ « Vp, r(,) 



(6) 



dt 



1 — C(i) 2e, 



M2 - PeE 



ret 
(i) 



(7) 



dsj 
'dt 



1 — C(i) 2e, 



dt 



r(i) dt 



(8) 



where the notation implies that the function is taken at the coordinates of the test- 
particle, i.e. = F{t,Xi{t),Pi{t),€i{t)). 

In Eqs.(6)-(8) the common multiplication factor (1 — C(j))~^ contains the energy deriva- 
tives of the retarded self energy 



C(i 



(0 



1 

2^ 



(9) 
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which yields a shift of the system time t to the 'eigentime' of particle i defined by ii — 
t/{l — C{i)). As in Refs. [41, 42] we will assume C{i) = furtheron. We mention, that 
the approximation C(i) = is also employed in the more standard transport models that 
operate with on-shcll quasi-particles. In the limiting case of particles with vanishing 
gradients of the width Txp the equations of motion (6) - (8) reduce to the well-known 
transport equations of the quasi-particle picture. 

Furthermore, following Ref. [40] we take = P'^ — ReT7^* as an independent variable 
instead of Pq, which then fixes the energy (for given P and M^) to 

^0 = + + ReTT^'p^, . (10) 
Eq. (8) then leads to 



dt r(j) dt 



11) 



for the time evolution of the test-particle i in the invariant mass squared as derived in 
Ref. [41]. It is worth noting that Eq. (11) (for r(j) 7^ 0) is equivalent to 

which states that the ratio of AM^ = Ml — Mq to the actual width r(j) is a constant in 
time. 

2.2 Collision terms 

The collision term of the Kadanoff-Baym equation in first order gradient expansion reads 
as [41] 

IcoiiiX, P, M^) = TT2TnTT^A{X, P, M^)A{X, P2, M^)A{X, P3, M|)A(X, P,, m|) 
\G{{P, M') + (P2, M^) ^ (Ps, Ml) + (P4, M^)) 1^ 5 5(4) (P + P2 - P3 - P4) 

[ ^XPsMl ^XPiMl fxPM^ fxP2Ml ~ ^XPM^ ^XP2M^ fxPsM^ f XP4,Ml \ i^"^) 

with 

fxPM'^ = ^+VNxPM^ (14) 

and f] = ±1 for bosons/fermions, respectively. The indices A,S stand for the antisym- 
metric/symmetric matrix element of the in-medium scattering amplitude G in case of 
fermions/bosons. In Eq. (13) the trace over particles 2,3,4 reads explicitly for fermions 



where cr2, T2 denote the spin and isospin of particle 2. In case of bosons we have 

since the spectral function is normalized as 

l'^A,(X,P) = l (17) 

whereas for fermions we get 

J^A,{X,P) = 1. (18) 

Here the spectral function Ap in case of fermions in Eq. (13) is obtained by considering 
only particles of positive energy and assuming the spectral function to be identical for 
spin 'up' and 'down' states. 

We note that the extension of Eq. (13) to inelastic scattering processes (e.g. —>■ 
Yn) or {nY — > KN etc.) is straightforward when exchanging the elastic transition 
amplitude G by the corresponding inelastic one and taking care of Pauli-blocking or 
Bose- enhancement for the particles in the final state. For bosons we will neglect a Bose- 
enhancement factor throughout this work since their actual phase-space density is small 
for the systems of interest. 

Thus the transport approach and the particle spectral functions are fully determined 
once the in-medium transition amplitudes G are known in their full off-shell dependence. 
For the transitions NN ^ A^A^ and NA <-> A^A^ we employ the same approximations as 
in Ref. [41] whereas for the strangeness sector KN ^ KN and KN <-> ttY we will use 
the off-shell transition rates as determined from the coupled channel G-matrix approach 
to be described below. 



3 The coupled-channel G-matrix approach 

In Ref. [29], the effective KN interaction in the nuclear medium (G-matrix) at tempera- 
ture T = was derived from a meson-baryon bare interaction built in the meson exchange 
framework [51]. As the bare interaction permits transitions from the KN channel to the 
ttE and ttA ones, all having strangeness S — —1, one is confronted with a coupled-channel 
problem. Working in an isospin coupled basis, the KN channel can have isospin 7 = 
or / = 1, so the resultant G-matrices are classified according to the value of isospin. For 
1 = 0, KN and vrS are the only channels available, while for 1 = 1 the ttA channel is also 
allowed. In a schematic notation, each G-matrix fulfills the coupled channel equation: 

(MiSi I G(0) I M2B2) = {MiBi I V{V^) I M2B2) 

+ Y: {MiB, I V{^s) I M^B^) J^^^^^ (M3E3 I G(Q) | M2B2) , (19) 
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where Q is the 'starting energy', given in the lab. frame, and y/s is the invariant center-of- 
mass energy. In Eq. (19), Mj and represent, respectively, the possible mesons {K, n) 
and baryons {N, A, E), and their corresponding quantum numbers, such as coupled spin 
and isospin, and linear momentum. The function Qm^b^ stands for the Pauli operator 
preventing the nucleons in the intermediate states from occupying already filled states. 

The prescription for the single-particle energies of all the mesons and baryons par- 
ticipating in the reaction and in the intermediate states is written - in nonrelativistic 
approximation but keeping relativistic kinematics - as 

EMi{B,){p) = ^Jp^ + mlj.^Bi) + ^M,(B,)(P, EZ.^Bi)iP)) > (20) 

where UMi{Bi) is the single-particle potential of each meson (baryon) calculated at the real 
quasi-particle energy E'^.^g.^. For baryons, this quasi-particle energy is given by 



Ef^{p) = ^p^ + ml + UBM , (21) 
while, for mesons, it is obtained by solving the implicit equation 

{E]^M"=P" + ^k + ^^^MME'^M ' (22) 

where S^* is the retarded meson self energy. 

The K single-particle potential in the Brueckner-Hartree-Fock approach (at tempera- 
ture T = 0) is given by 

Uk{Pk, ) = E i^N I Gkn-.kn{^ = E'S + £;f ) I KN) , (23) 

N<F 

where the summation over nucleon states is limited to the occupied nucleon Fermi sphere. 
The K self energy is obtained from the optical potential (neglecting energy derivatives) 
through the relation 



WiPR, ^) = ^\/Pr + < Uk{Pk, ^) , (24) 

with uj = E'^^{pj^). As it can be easily seen from Eq. (23), since the KN effective 
interaction (G-matrix) depends on the K single-particle energy, which in turn depends 
on the K potential through Eqs. (22), (24), one needs to solve a self-consistent problem. 

After self-consistency is reached, the complete energy-and-momentum dependent self 
energy of the K can be obtained from Eq. (24), which allows one to derive the K propa- 
gator 

and the corresponding spectral density, defined as 

Sk{Pk,^) = --^'^E)k{pk,u) . (26) 

TT 
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We note that our self-consistent procedure amounts to replace the complex energy de- 
pendent self energy, S^*(p^,a;), in the K propagator by that evaluated at the quasi- 
particle energy, = E'^{pj^)). We denote this procedure as the quasi-particle 

self- consistent approach, which retains the position and the width of the peak of the K 
spectral function at each iteration. This simplification relative to using the complete en- 
ergy dependence, as done in Refs. [27, 28], allows one to perform analytically the energy 
integral of the intermediate loops, thus reducing the four- dimensional integral equation 
to a three-dimensional one. 

An important modification for the antikaon properties comes from considering the 
pion self energy, H'^^^^ {p^^ , uj) , in the intermediate vrS, ttA states present in the construction 
of the KN effective interaction. Since the dressing of the pions has a strong influence 
on the antikaon in-medium properties (cf. Refs. [28, 38] and also Sections 3.2 and 5 of 
the present work), it is of relevance to be more specific on the actual realization. The 
pion self energy of the present work is built up from ph, Ah and 2p2h excitations and 
contains the effect of nucleon-nucleon short-range correlations via a phenomenological 
model developed in Refs. [52, 53], which gives rise to a smooth momentum-dependent 
Landau-Migdal parameter g' of the order of 0.6. This is very similar to the value obtained 
from microscopic calculations based on the nucleon-nucleon G- matrix [54, 55, 56] including 
the crossed-channel contributions, similarly as the induced interaction of Babu and Brown 
[57]. Moreover, the phenomenological pion self energy used here has been tested in various 
reactions involving the interaction of virtual or real pions with nuclei in Refs. [58, 59, 60, 
61]. Certainly, this pion self energy is extrapolated to higher densities in the present work 
and it might not be justified to use the same value for g'. However, the microscopic study 
of Ref . [54] , performed in the context of the problem of pion condensation at high baryon 
density, only showed a weak dependence of the nucleon-nucleon residual interaction with 
density. Nevertheless, we regard the issue of pion dressing at higher baryon densities still 
to be open and will present our G-matrix results in two limits, i.e. with and without pion 
dressing. 

In order to address energetic nucleus-nucleus collisions, the coupled G-matrix equa- 
tions in principle have to be extended to non-equilibrium phase-space configurations. This 
is quite a formidable task such that we restrict to finite temperature calculations since 
the antikaon dynamics of interest proceeds in the environment of a hot hadronic fireball. 

The introduction of temperature in the G-matrix equations affects the Pauli blocking 
of the intermediate nucleon states as well as the dressing of mesons and baryons. The 
G-matrix equation at finite T reads formally as in Eq. (19), but replacing (cf. Ref. [38]) 

Qmb Qmb{T) 
G{Q) G{Vl,T) 
Em , Eb — > Em{T) , Eb{T) . 

The function Qmb{T) is unity for meson- hyperon states while, for KN states, it follows 
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Qmb{T) = 1 — n{pN,T,ij,) with the nucleon occupation number given by 



1 + exp 



T 



(27) 



The chemical potential /i in Eq. (27) is obtained by imposing the normalization condition 
for the nucleon density 

P = j d^PN n{pN, T, n) , (28) 

where = 4 is the degeneracy factor for symmetric nuclear matter. 

A finite temperature also affects the properties of the particles involved in the process. 
The K optical potential at a given temperature T is calculated according to 

URiPR, Ef, T) = j d'pN nipN, T) {KN \ Gr^^r^{SI = E'j^ + , T) | KN) , (29) 

which again is a self- consistent problem for Ur. More explicitly, using the partial wave 
components of the G-matrix, we obtain 

Uk{pk,E'^IT) =1 ^(2J+1)(2/+1) [ nipN,T) p% dpN (30) 

^ L,J,I 

where p and P are the relative and center-of-mass momentum, respectively, averaged over 
the angle between the external K momentum in the lab system, pR, and the internal 
momentum of the nucleon, pjv- In the actual calculations, we include partial waves up to 
L = 4. In extension of the work in Ref. [38], we here additionally include the S*(1385) 
resonance dynamics in a similar way as in Ref. [23]. We note that the role of this sub- 
threshold E*(1385) resonance is neghgible for KN free scattering observables but, as 
suggested in Ref. [39], it can become relevant in the medium since, due to the attraction 
felt by the antikaons, one is effectively exploring lower values of y/s. 



3.1 Antikaon quasi-particle properties 

With respect to nucleus-nucleus collisions at SIS energies (1-2 A-GeV), we have performed 
the G-matrix calculations at a fixed temperature T = 70 MeV which corresponds to an 
average temperature of the 'fireballs' produced in these collisions. Moreover, we note that 
variations in the temperature from 50 - 100 MeV do not sensibly affect the quasi-particle 
properties in the medium [38]. 

Whereas in Ref. [38] the antikaon potentials have been shown at the quasi-particle 
energy (as a function of momentum and nuclear density), we display in Fig. 1 the real 
part of the antikaon potential Re Ur Si function of 

V~s = - Ph (31) 
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where uj denotes the antikaon energy and its momentum relative to the nuclear mat- 
ter rest frame, for different densities and momenta Px= 0, 150, 300, and 500 MeV/c, 
respectively. The arrows in Fig. 1 show the pole mass of the antikaon in free space for 
orientation. It is seen from Fig. 1 that the real part of the potential is attractive through- 
out and approximately linear in the nuclear density from 0.5 po to 2 p^. Furthermore, the 
potential is slightly more attractive for smaller off-shell masses than for masses above the 
free invariant mass. With increasing momentum the potential becomes shallower and 
the dip at an invariant mass of ~ 320 MeV vanishes. We note, that the results presented 
in Fig. 1 stem from a G-matrix calculation including the pion dressing. When excluding 
pion self energies the antikaon potential becomes even more attractive especially for low 
invariant masses. 

The imaginary part of the antikaon potential Im is shown in Fig. 2 also as a 
function of ^ for different densities and momenta = 0, 150, 300, and 500 MeV/c, 
respectively. The imaginary part is also roughly proportional to the nuclear density and 
larger for smaller invariant masses than for masses higher than the free antikaon mass 
(arrows). The structure observed around ~ 300 MeV at low momentum is due to 
resonant p-wave coupling to SiV^^ states as will be discussed below. In general, the 
imaginary part shows only a weak dependence on the momentum . 

The resulting spectral function Sk{px,u!) is shown in Fig. 3 as 



where ^/s can be identified with the invariant mass M, for different densities and momenta 
Px = 0, 150, 300, and 500 MeV/c, respectively. In line with Figs. 1 and 2 the maximum of 
the spectral function shifts to lower invariant masses with increasing density and becomes 
substantially broader. The relative change with momentum is only moderate as seen 
from Fig. 3. 

As pointed out in Section 2, the information presented in Figs. 1 and 2 - on a much 
finer grid in momentum and density - is used for the off-shell propagation of antikaons in 
the nuclear medium according to the Eqs. (6), (7), (8) in the transport approach. On the 
other hand, the spectral functions from Fig. 3 - also on a much finer grid in momentum 
and density - enter the collision terms of Eq. (13) for the production of antikaons in 
nucleon-baryon, meson-baryon or pion-hyperon interactions. 

3.2 In-medium transition rates 

Apart from the spectral functions Axp the collision terms in Eq. (13) depend on the local 
phase-space densities N^pj^2-i that are calculated dynamically for all hadrons, and the 
transition matrix-elements squared |G'((P, M^) (-P2,M|) (-P3,M|) + {PA,Ml))\\g. 
Here the latter have to be known not only for on-shell but also for off-shell particles. 

For binary reactions involving antikaons like KN — > KN or ttA these matrix 

elements are determined by the C-matrix equation (19) such that no new parameters 




(32) 
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or unknown cross sections enter the transport calculations. Actual cross sections are 
determined as a function of the invariant energy squared s as 

.„,_3..(.) = (2.)=£iM5£i^ _J__^^ a^a, (33) 

where p and p' denote the center-of-mass momentum of the particles in the initial and final 
state, respectively. The sums over i and a indicate the summation over initial and final 
spins, while si,S2 are the spins of the particles in the entrance channel. Apart from the 
kinematical factors, the transition rates are determined by the angle integrated average 
transition probabilities 

= /dcosW ^^_J^^X:X; GtG (34) 

which are uniquely determined by the G-matrix elements evaluated for finite density p, 
temperature T and relative momentum with respect to the nuclear matter rest frame. 

Before coming to the actual in-medium problem we show in Fig. 4 a comparison of 
our calculations for the strangeness exchange cross sections K~p — > E°7r°, K~p E+vr", 
K~p Att^, and K~p S^7r+ in 'free space' with the corresponding experimental data 
from [62] as a function of the antikaon momentum p^ in the laboratory. As seen from 
Fig. 4 we miss the small resonance for ^ 0.4 GeV/c in the K~p — > S°7r° channel^ and 
slightly overestimate the K~p — > Att^ channel. Otherwise, the comparison shows that our 
transition probabilities arc sufficiently realistic for the vacuum cross sections. It should 
be pointed out, that the A(1405) is generated dynamically in our approach whereas in 
Refs. [63, 64] it is treated as an elementary field. 

We now turn to the transition probabilities of Eq. (34) at finite density p, finite 
temperature T and finite antikaon momentum in the nuclear matter rest frame. They 
are displayed in Figs. 5, 6, 7 and 8 for the reactions K~p K~p, K~p S°7r°, K~p 
Att^ and Att^ Avr", respectively. The isometric plots show the probabilities of Eq. (34) 
as a function of density p (in units of po) and ^/s for a momentum = at T = 70 
MeV. The results on the l.h.s. correspond to a G-matrix calculation without pion dressing 
whereas those on the r.h.s. include pion dressing as described in Section 3. We note that 
the general shape of these transition probabilities does not change very much at finite 
momentum pf^ such that we discard an explicit representation over this variable. 

In order to pick up the physics from Figs. 5-8 we recall that the K~p threshold in free 
space corresponds to ~ 1.432 GeV, whereas the thresholds for 7r°A and 7r°E are at 
1.25 GeV and 1.332 GeV, respectively. In free space (p = 0) the coupling to the A(1405) 
resonance provides the dominant matrix elements below (and close to) the KN threshold 
in all channel amplitudes that contain the / = component (Figs. 5 and 6). Note that, 
in the case of off-shell antikaon dynamics in the nuclear medium, also lower invariant ^/s 
become accessible such that a resonant amplitude slightly below the free threshold will 
give large in-medium transition rates. The same considerations apply to the p-wave E* 

^The A(1520) is not included in our calculations. 
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resonance located at 1.385 GeV. Its influence is seen more clearly in the 1=1 channel 
Att^ — > Att" (Fig. 8). The reason is that, since the p-wave amplitude is proportional to the 
product of center-of-mass momenta in the incoming and outgoing channels, pp', it gets 
enhanced for invariant energies above and away of the corresponding thresholds. This also 
explains the enhanced contribution of this resonance in the channel K'p Air^ (Fig. 7) as 
density grows due to the fact that the K~p threshold, which lies above the S*(1385) in free 
space, also moves to lower energies in the medium. This is especially visible on the l.h.s of 
Fig. 7 when pion dressing is ignored. In fact, the width of the E*(1385) resonance already 
increases in the medium when only the dressing of the antikaons is incorporated due to 
the opening of new decay channels, such as H*N AN, UN, nAN, nUN. However, the 
additional dressing of the pions in the vrA, vrE G-matrix intermediate states, to which the 
E* couples very strongly, enhances tremendously its decay width into the above mentioned 
nucleon-induced channels, hence producing a very 'smeared out' contribution, as observed 
on the r.h.s. of Fig. 7. 

We also note that in the region of the vrS threshold the squared matrix amplitudes 
that contain / = components (see Figs. 5 and 6) present a double peak structure at 
low densities when pion dressing is included. This was already pointed out in Ref. [38], 
where it was shown that the in-medium 7 = s-wave resonance moves down below the 
ttE threshold at normal nuclear matter density po- Hence, as density increases and, for 
some range of densities, this resonance shows up as two distinct bumps. 

In any case, it is clear that, whatever resonant structures are present at low densities, 
they melt away already at a modest density of ~ 0.5po for all calculations that include 
pion dressing (r.h.s.). 

Finally, one observes very modest structures, especially visible in the channels that 
contain 1 = 1 contributions at low energies (sec Figs. 7 and 8), which explain the behaviour 
found in the antikaon potential (cf. Figs. 1 and 2) as well as the bump of the antikaon 
spectral function (cf. Fig. 3) around ^/s ~ 300 MeV, i.e. roughly 200 MeV below the 
free K mass. These structures correspond to the couphng of the K meson to E-nucleon 
hole configurations. 

4 Antikaon production from NN and tvN collisions 

Whereas the G-matrix approach described in the previous subsection uniquely determines 
the off-shell transitions KN — > KN and KN ttY, i.e. the dynamics of hadrons with 
s-quarks in the medium, the production cross sections of kaons and antikaons from NN 
or ttN coUisions are not specified accordingly. As already shown in Ref. [65] the latter 
channels are subdominant relative to the leading s-quark exchange reactions ttY — > KN. 
This dominance of the nY reaction channel in nucleus-nucleus collisions at SIS energies 
is related to the fact that kaon cross sections and (due to strangeness conservation) the 
hyperon cross sections arc about 2 orders of magnitude higher than the antikaon cross 
sections [65, 66]. Consequently ss-quarks are essentially created with KY final channels 
from NN, AN and nN reactions and a large fraction of antikaons stems from the strange 
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flavor exchange channels ttY — > KN. We mention that with increasing bombarding energy 
also the channel rcn KK may contribute significantly for heavy systems [67, 68]. 

For the determination of the in-medium NN — > NNKK and 'kN — > NKK cross 
section we here proceed as follows: from experimental data on the 'free' production cross 
sections we extract an average matrix element squared \M.\^ by dividing out the phase 
space for on-shell particles in the final state and the fiux factor (cf. Ref. [69]). Using 
the same matrix element, the off-shell cross sections then are obtained by employing 
the dynamical spectral functions for the hadrons and correcting for the modified final 
state phase-space for off-shell hadrons. The actual implementation is done as follows: In 
case of meson production by off-shell baryon-baryon or meson-baryon collisions we either 
have 2 (e.g. vrA^ ^ K+K/T.), 3 (e.g. NN K+AN or K+J:N) or 4 particles (e.g. 
NN — > NNK^ K~) in the final channel. Since the final mesons may be off-shell as well, 
one has to specify the corresponding mass-differential cross sections that depend on the 
entrance channel and especially on the available energy ^/s in the entrance channel. 

We start with the explicit paramctrizations for meson (m) production cross sections 
given in Ref. [7] for on-shell mesons as a function of the invariant energy ^/s in case of 
nucleon-nucleon or pion-nucleon collisions, i.e. CNN^NNmiVs) or (TnN^mN{y/s), respec- 
tively, that are sufficiently well controlled by experimental data in 'free' space. To this 
aim we show in Fig. 9 the respective data for pp — > K'^K^pp, pp — > K^K~pp from Refs. 
[62, 70, 71] and 7r~p K~K^p from Ref. [62] in comparison to our paramctrizations. The 
reader should not worry about the fact that the phase-space oriented parametrization for 
the pp reactions does not match the 2 lowest points close to threshold since the latter are 
enhanced by the strong final state interaction in the pp and Kp final channel in free space. 
On the other hand, such final state interactions are essentially 'screened' in the nuclear 
medium (especially at high density) such that the phase-space oriented approximations 
should hold sufficiently well. 

Keeping this concept in mind, the in-medium mass differential cross sections - far 
above the corresponding thresholds - are approximated by 

— = (JNN-^mNNiJs - J ) Am{M ,P,p), (35) 



where AmiM"^ ,p. p) denotes the meson spectral function for given invariant mass M^, 
relative momentum p and nuclear density p as determined from G-matrix theory in Section 
3. In Eq. (35) the threshold energy — Mq + Ml -\- depends on the masses of the 
hadrons in the final channel, i.e. Mo,Mj* and M|. Actual events then are selected by 
Monte-Carlo according to Eq. (35). Close to threshold ^/sj, i.e. for a/s — Mq — M^* — M2 < 
2Ttot-, where M^*, denote the final off-shell masses of two nucleons, Mq the meson pole 
mass and Vtot its total width, the differential production cross section is approximated by 
a constant matrix element squared times the available phase-space, 

da^^^rnM^s) ^ \Mrr.?A^{M\p,p) R,{s, , M^,* , M^) , (36) 

where the matrix element \A4rn\ is fitted to the on-shell cross section typically from 50 to 
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500 MeV above threshold. In Eq. (36) the function denotes the 3-body phase-space 
integral [72] in case of a mNN final state. 

In case of 4 particles in the final state, i.e. in the channel NN — > KKNIN2, where 
the A^*'s denote off-shell nucleons, the differential cross section is approximated by 

E E E E ^^^^-BB^ Af ATMi Ma ( V^) ^ /g^-j 

d^P\ d?p2 d^Ps d^p4 
ctbb^nnmi M2{\/s) ^„„^. (Pi + P2-P1-P2-P3 -P4), 

where R4 denotes the 4-body phase-space integral [72]. Similar strategies have been 
exploited in case of subthreshold KK or even pp production in proton-nucleus and nucleus- 
nucleus collisions in Rcfs. [41, 65, 73, 74]. 

The resulting cross sections for K~ production from NN and 7r~p collisions using 
on-shell nucleons and kaons in the final state, however, employing the antikaon spectral 
functions from the G-matrix approach in Section 3, are displayed in Fig. 10 as a function 
of ^/s for different nuclear densities ranging from 0.25 po to 2.25 po in comparison to the 
cross section in free space (solid lines). With increasing nuclear density the subthreshold 
production of mesons becomes enhanced considerably relative to the respective vacuum 
cross section due to a shift of the antikaon pole mass and a broadening of its spectral 
function (cf. Fig. 3), but the absolute magnitude stays small below threshold even for 
p ~ 2po- We recall again that the ttN and NN production channels for KK pairs arc not 
the leading channels at SIS energies and minor uncertainties in the off-shell treatment of 
the production cross sections are unlikely to show up in the final K~ abundancies and 
spectra [65, 66]. 

We use isospin symmetry to relate the cross sections for Tr~p or pp induced reactions to 
ttN and A^A^ collisions. Furthermore, due to a lack of any experimental information, the 
ttA production channel is assumed to be same as the ttN channel at the same invariant 
energy y^. The NA and AA production channels again are taken the same as the NN 
channel for fixed invariant energy. Apart from interactions between baryons and pions 
with baryons also meson-meson collisions become important with increasing bombarding 
energy [67, 68]. We include the tttt ^ KK channel using the cross section from Ref. [7] 
and in addition to Ref. [65] also the channel n + K ^ K* employing a Breit-Wigner 
resonance cross section with K* resonance properties from Ref. [75]. Moreover, the 
decay —>■ KK is included in the transport calculation which, however, is found to be 
subdominant. 

Now all matrix elements, self energies and/or cross sections for antikaon production 
and propagation are specified such that we can continue with the actual transport cal- 
culations. We recall, that in case of kaons or antikaons at SIS energies (or charmonia 
at SPS energies) we treat the latter hadrons perturbatively as in Ref. [65, 76, 77], i.e. 
each testparticle achieves a weight Wi defined by the ratio of the individual production 
cross section to the total ttB or BB cross section at the same invariant energy. Their 
propagation and interactions are evaluated as for baryons and pions, however, the baryons 
(pions) are not changed in their final state when interacting with a 'perturbative' particle. 
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In case of strange flavor exchange reactions the individual weight of a -ft^ is given to the 
hyperon Y and viceversa. In this way exact strangeness conservation can be achieved 
during the transport calculation while obtaining reasonable statistics also for antikaons 
at 'subthreshold' energies. We mention again that the actual (space-time dependent) 
antikaon width is determined by the G-matrix approach described in Section 3. 

5 Nucleus-nucleus collisions 

We carry out the concrete applications for nuclear reactions at SIS energies (1.5 - 2 A 
GeV) that have been analysed within conventional transport models to a large extent (cf. 
Refs. [7, 16] and Refs. therein). However, before coming to the actual results we hke 
to address the off-shell propagation of antikaons and the evolution of the K~ spectral 
density in time and invariant mass for central nucleus-nucleus collisions. To this aim we 
show in Fig. 11 the time evolution of the quantity 

F{M, t)^ j dh j d^p A(r, t, p) A^(r, t, M^ p) (38) 

as given by the off-shell test-particle approximation [Eq. (5)]. The isometric plots in Fig. 
11 show F{M,t) for central collisions of Ni + Ni at 1.8 A-GeV for antikaons stemming 
from baryon-baryon (BB) (r.h.s.) and nY reactions (l.h.s.), separately. It is seen from 
Fig. 11 that the initial distribution of antikaons is widely spread in invariant mass for 
both production channels and that - with the expanding nuclear system - the antikaons 
become practically on-shell for large times, i.e. when moving to 'free' space with the 
expanding fireball. In case of the BB production channel the initial spreading in mass 
also shows sizeable contributions for M > 0.5 GeV whereas for the nY flavor exchange 
reaction dominantly low mass antikaons are produced which are kinematically not allowed 
in case of 'free' antikaon masses. The energy - to get on-shell flnally - stems from the 
collective expansion of the hadronic system. 

5.1 K"^ spectra at SIS energies 

Since kaons couple only weakly to nucleons and are not absorbed at low energies their 
collisional width is rather small such that they may be treated on-shell to a good approx- 
imation as in Ref. [77]. We recall, that the production channel NA NK~^Y, where 
Y denotes a hyperon, as well as the AA — > K^NY channel is not known experimentally 
and simple isospin factors as extracted from pion exchange [77] might not be appropri- 
ate. Though there are some efforts to resolve this uncertainty within extended boson 
exchange models [78], the latter models will hardly be tested experimentally. We note 
that the iVA — > NK^Y and AA — > K^NY cross sections employed in the transport 
models of Refs. [79, 80] are larger by factors of 2 - 3 than ours, which leads to higher 
K'^ (and hyperon) cross sections from nucleus- nucleus collisions in the 'free' scenario [7]. 
These cross sections are reduced again in Refs. [79, 80] by a repulsive kaon potential in 
order to achieve a better agreement with measured kaon cross sections. 
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In this work we do not address this particular question in more detail since only the 
kaon abundancies are of interest here due to the associated productions with hyperons Y — 
A, E. As mentioned above, the channel nY — > K'N is expected to be dominant such that 
the transport calculations have to reproduce the (experimental) hyperon abundancy with 
sufficient accuracy. Since there arc no explicit hyperon spectra available from nucleus- 
nucleus collisions at SIS energies, we will perform a detailed comparison to measured 
spectra (see below). 

On the other hand, antikaons couple strongly to nucleons and thus achieve a large 
coUisional width in the nuclear medium as demonstrated in Fig. 2. Accordingly, off- 
shell antikaons might be produced at far subthreshold energies (cf. Fig. 10), become 
asymptotically on-shell (cf. Fig. 11) and thus enhance the K~ yield. In fact, as shown 
in the model study in Ref. [41] the K~ yield might be enhanced up to a factor of 2 for 
Ni + Ni at 1.8 A-GeV when including the antikaon off-shell propagation. However, in 
the latter study the off-shell transition rates nY K~N had been extrapolated from the 
on-shell rates in free space which according to Figs. 5-8 should be questionable. 

In the following we will show three different limits simultaneously in comparison to 
the experimental data: i) a calculation without any antikaon in-medium effects (denoted 
by 'free'), where the G-matrix elements at density p — are adopted and the antikaon 
spectral function is taken as a 5-function on-mass shell; ii) a calculation using the G- 
matrix elements, spectral functions and potentials from a G-matrix calculation without 
pion dressing and iii) a calculation using the G-matrix elements, spectral functions and 
potentials from a G-matrix calculation with pion dressing. We mention, that we use only 
the G-matrix elements calculated at an average temperature T— 70 MeV and do not 
follow the change of these matrix elements with decreasing temperature T. Apart from 
the tremendous numerical effort to calculate the G-matrix additionally on a narrow grid 
in temperature T, the modifications with temperature are rather moderate or even small 
(cf. Ref. [38]), such that we discard this variation. 

The actual spectra for the systems G + G at 1.8 A-GeV and 9cm = (90 ± 10)", 
Ni+Ni at 1.8 A-GeV and 9iat> = (44±4)° and Au+Au at 1.5 A-GeV and 9cm = (OOilO)" 
are shown in Figs. 12-14 in comparison to the data from Refs. [81, 82, 83, 84]. The K'^ 
spectra are shown in the upper plots whereas the K~ spectra are displayed in the lower 
parts, respectively. In Figs. 12-14 the spectra are displayed for the full G-matrix 
calculations only, since the spectra are practically insensitive to the different limits 
addressed for the antikaons within the statistics of the transport calculations. 

In case of G -|- G (Fig. 12) the 'slope' Tq of the K'^ spectrum, defined by 

is shghtly underestimated when comparing the data at 9cm = (90 ± 10)° with the calcu- 
lations in the same angular bin (solid line). However, when averaging the calculated K'^ 
spectrum over the solid angle Q, in the center-of-mass system (dashed line) the experimen- 
tal spectrum is well described in shape as well as absolute magnitude. This comparison 
demonstrates that the details of the double differential cross with respect to the kaon 
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kinetic energy in the cms, E^^^ and angle Vt are not fully reproduced by the transport 
calculation, but the kaon abundancy and average spectrum compare reasonably well. Due 
to the associated production of hyperons (with kaons) we thus conclude that also the av- 
erage hyperon abundancy and spectra should compare sufficiently well with experiment 
though explicit data arc not available for the A and S spectra. 

For Ni + Ni at 1.8 A-GeV (Fig. 13) the slope of the spectrum as well as the 
abundancy is slightly too low for the calculations at Oiab = (44 ±4)° (solid line) in 
comparison to the data in the same angular range. Again a very satisfactory description 
is obtained in case of the angle averaged calculated spectrum (dashed line) as in case of 
the C + C system in Fig. 12, such that the interpretation and conclusions drawn in the 
context of the latter system also hold for Ni + Ni at 1.8 A-GeV. Furthermore, for Au + Au 
at 1.5 A-GeV (Fig. 14) the multiplicity seems to be roughly in line for the spectra in 
the angular range 9cm — (90 ± 10)° (sohd hne) as well as for the angular averaged spectra 
(dashed line) in comparison to the preliminary data from Ref . [84] . 

The KaoS [85] and FOPI Collaborations [86] have independently measured pro- 
duction in Ni + Ni collisions at 1.93 A-GcV. In fact, their experimental results for the 
inclusive rapidity distributions agree quite well in the common region of acceptance in 
rapidity such that systematic experimental uncertainties are much better under control. 
In Fig. 15 we compare our calculated rapidity distributions for this system with the data 
from Ref. [85] for 'semi-central' collisions that correspond to impact parameter h < 4.5 
fm. For Ni + Ni at 1.93 A-GeV our calculations almost perfectly agree with the mea- 
sured rapidity distribution since deviations in the double differential spectra with 
respect to momentum and angle do not show up any more in the rapidity distributions. 

The K'^ spectra from our transport calculations thus are found to be in 'approximate' 
agreement with the data for the various systems, however, systematic uncertainties in the 
order of 30% in the (and associated hyperon abundancies) cannot be excluded. On 
the other hand, it is presently not clear if the different data sets are compatible with each 
other to a higher accuracy. 

We now turn to the results for the K~ spectra (lower parts in Figs. 12-15). Whereas 
for the 'free' calculations (dashed lines) the experimental K~ spectra are substantially 
underestimated for C + C and Ni + Ni at 1.8 A-GeV (in line with the calculations in Ref. 
[7, 65]) the preliminary data for Au + Au at 1.5 A-GeV in Fig. 14 appear to be described 
rather well in magnitude. Only the slope of the calculated K~ spectrum is slightly too 
hard in this limit. Surprisingly, the K" rapidity distribution for semi-central collisions 
of Ni + Ni at 1.93 A-GeV (Fig. 15) is underestimated only by less than a factor of 2 in 
the 'free' case. This result - on first sight - appears incompatible with the comparison 
displayed in Fig. 13 for Ni + Ni at 1.8 A-GeV where the K~ spectra are underestimated 
on average by more than a factor of 4 in the 'free' case^. 

The calculated slope Tq of the K~ spectrum is too high for the C + C system, roughly 
comi)atil)le for Ni + Ni at 1.8 A-GcV and slightly too high again for Au+Au at 1.5 A-GeV 

^According to private communication with the KaoS Collaboration the normalization of the K~ data 
for Ni + Ni at 1.8 A-GeV is estimated to be too high by up to a factor of 2, which is on the lower level 
of the error bars quoted in Ref. [82] . 
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in the 'free' scenario when compared to the data in Figs. 12-14. The G-matrix calculations 
including pion dressing (full lines with open triangles) for C + C at 9cm = 90° ± 10° are 
approximately compatible in magnitude with the 'free' calculations for kinetic energies 
above 100 MeV, however, the slope of the spectrum is now in better agreement with the 
data. For Ni + Ni at 1.8 A-GeV and 6iab = 44° ± 4° the G-matrix calculations even 
fall slightly below the 'free' result whereas for Au + Au at 1.5 A-GeV they are again 
compatible with the data and the 'free' result within statistics for > 100 MeV. This 
also holds for the K~ rapidity distribution for Ni + Ni at 1.93 A-GeV in Fig. 15 where 
the full G-matrix calculations reproduce the experimental spectrum almost 'perfectly'. 

In general, the medium modifications of the K~ spectra arc found to be only mi- 
nor within the G-matrix calculations including pion dressing in the regime of antikaon 
momenta Pcm or kinetic energies E^^, where explicit data are available. This is sur- 
prising since the antikaon spectral functions show a sizeable shift of strength to lower 
invariant masses (Fig. 3) such that - according to phase-space - the cross sections should 
be enhanced. However, this expectation is only valid if the transition and production 
amplitudes do not change in the medium, which actually does not hold for the present 
G-matrix calculations as seen from Figs. 5-8. Accordingly, the dominant production 
channel n + Y ^ KN decreases strongly with nuclear density due to the rapid melting 
of the A(1405) and E(1385) resonances with density such that less s-quarks now can be 
transferred from hyperons to antikaons at high baryon density. Qualitatively, our findings 
are similar to the self-consistent results of Ref. [37] . 

The relative role of the hyperon 'resonances' becomes more clear when looking at 
the results for the G-matrix calculations without pion dressing (solid lines with open 
circles in Figs. 12-15). In this case all spectra are found to be enhanced relative to the 
'free' calculations and the G— matrix calculations with pion dressing. The experimental 
spectra for G + G at 6cm = 90° ± 10° are still underestimated as well as for Ni + Ni at 
1.8 A-GeV and Oiab = 44° ± 4°. However, the calculations for Au + Au at 1.5 A-GeV 
and Ni + Ni at 1.93 A-GeV now are clearly above the measured data. In hne with 
the higher K~ multiphcity also the spectral slope softens shghtly because the antikaon 
spectral strength is shifted to lower invariant masses in case of the 'surviving' subthreshold 
resonances A(1405) and S(1385), respectively. We argue, that the calculations with and 
without pion dressing represent a 'band of uncertainty' within the G-matrix calculations, 
that implicitly depend on the interaction schemes involved. Nevertheless, in all limits 
considered here there is no convincing description of all spectra simultaneously! 

The question thus arises if medium effects for antikaons show up in kinematical regimes 
not accessible to the present detector setups. To this aim we show in Fig. 16 the Lorentz 
invariant K' spectra ior C + C and Ni + Ni at 1.8 A-GeV and Au + Au at 1.5 A-GeV 
as a function of the antikaon momentum Pcm in the cms frame. Whereas above ~ 0.5 
GeV/c all spectra within the hmits addressed before are roughly comparable, a large 
enhancement at low momenta is found for all systems within the G-matrix calculations 
relative to the 'free' case. This enhancement is even more pronounced when excluding 
pion dressing (full lines with open circles). The K~ enhancement at low momenta is 
essentially due to a deceleration of antikaons in the combined attractive Coulomb and 
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nuclear mean field as well as due to the off-shell propagation of 'low mass' antikaons. In 
the off-shell propagation (see Eqs. (6)-(8)) spectral components with masses below the 
free pole mass (cf. Fig. 11) decrease in momentum during the expansion phase according 
to Eq. (7) in order to become on-shell in free space. 

The results from Fig. 16, which qualitatively agree with the earlier analysis by Wang 
et al. [87] using on-shell transport calculations, also explain to some extent why the cal- 
culations within the full G-matrix approach are compatible with the K'^ rapidity spectra 
for Ni -\- Ni at 1.93 A-GeV (Fig. 15), that experimentally have been obtained by extra- 
polation to low antikaon momenta, but not with the K~ spectra for Ni + Ni at 1.8 A-GeV 
(Fig. 13) in the higher momentum range of the KaoS acceptance. The reason is most 
likely related to the strong increase of the antikaon spectrum at low momenta in the cms. 

5.2 / ratios versus centrality 

We continue with the / ratio in Au + Au collisions at 1.5 A-GcV as a function of 
centrality which has been also addressed experimentally by the KaoS Collaboration [84]. 
In Fig. 17 we show this ratio as a function of the number of participating nucleons Apart 
for the three limits discussed before including a cut Pcm > 0.35 GeV/c which roughly 
corresponds to the acceptance of the KaoS spectrometer in Ref. [84]. In the 'free' calcula- 
tions as well as G-matrix calculations without pion dressing this ratio - within statistics 
- increases with centrality or Apart- This increase is no longer present for the G-matrix 
calculations with pion dressing (full line with open triangles). The preliminary data for 
this ratio from the KaoS Collaboration [84] (full circles) show - within errorbars - a shght 
decrease with Apart or an approximately constant value of ~ 2%. This experimental value 
is not compatible with the G-matrix calculation excluding pion dressing as well as the 
'free' calculation, however, in the order of the full G-matrix results. 

Some comments on the / ratio appear necessary since the dependence on cen- 
trality and its actual values strongly correlate with the angle and momentum cut intro- 
duced. This comes about because the K~ spectra show a softer slope in the momentum 
spectrum relative to the mesons, which wc address to attractive in-medium effects 
for the antikaons. An alternative or complementary interpretation persists in relating the 
different slopes to different 'freeze-out' times since mesons decouple early in view of 
their low cross section with baryons, whereas antikaons interact strongly up to rather low 
densities of the expanding hadronic fireball. Consequently the relative change of the K" 
to spectral slopes with centrality strongly affect the K~ / ratio in Fig. 17 (for a 
high momentum cut) such that no 'simple' interpretation can be drawn directly. Without 
explicit representation we mention that in our calculations the K~ /t:~ and K'^ ji:^ ratios 
increase both with centrality (or Apart)- However, also the difference in the slopes for 
and K~ mesons increases with Apart for the full G-matrix calculations including pion 
dressing. 
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5.3 angular distributions 

We finally come to the K"^ angular distributions in the cms for Au + Au collisions at 1.5 
A-GeV that also have been measured by the KaoS Collaboration [88] for two different cen- 
trality cuts. In Fig. 18 we show these angular distributions for the three limits discussed 
before for K" mesons (lower part), whereas the angular distributions for K'^ mesons (up- 
per part) are shown for the full G-matrix calculations, only, since the same distribution 
has been obtained in the three limits (within statistics). All angular distributions are 
normalized to unity for cos^^cm = 0. 

We see that the angular distributions from our calculations for semi-central (6 < 
6 fm) and non-central (6 > 6 fm) collisions are more isotropic than the data. This points 
towards a lower amount of K'^ rescattering in the data [88] especially for non-central 
collisions. On the other hand, all three limits are compatible with the experimental 
measurements for i^~-mesons in semi-central reactions, which are roughly compatible 
with an isotropic distribution within statistics. Differences appear only for non-central 
reactions, where the 'free' (dashed line) and G-matrix calculation without pion dressing 
(sohd line with open circles) are closer to isotropy while the full G-matrix calculations 
with pion dressing (solid lines with open triangles) show a distribution, that is forward- 
backward peaked in the cms in better agreement with the data. We attribute this result 
to the lower amount of s-quark exchange reactions in the full G-matrix calculations (cf. 
Section 3). 

6 Summary 

In this work we have studied the production and propagation of antikaons with dynamical 
spectral functions Ah{X, P, M^) using off-shell transport theory. The in- medium proper- 
ties of the antikaons have been determined in a coupled-channel G-matrix approach in- 
cluding elastic scattering, charge exchange and s-quark exchange reactions with baryons. 
Antikaon mean-field potentials as well as spectral properties are uniquely determined 
within the G-matrix approach. However, the actual results strongly depend on the many- 
body scheme involved, especially on the dressing of the pions. Though in all cases the 
spectral strength of the antikaon is shifted to lower invariant masses with increasing den- 
sity, we find that, when including pion self energies, the A(1405) and S(1385) resonance 
structures in the transition probabilities melt away with baryon density already at ~ 0.25 
Po- This implies that K~ absorption as well as production from pion-hyperon collisions 
should be strongly suppressed in nuclear media as produced in nucleus-nucleus collisions 
at SIS energies (1-2 A-GeV) or even proton- nucleus reactions at 2-3 GeV [89]. 

From our dynamical calculations we find that the experimental spectra for ^^G 
C and ^^Ni -\-^^ Ni are underestimated in the 'free' as well as full G-matrix transport 
approach whereas the preliminary spectra for Au + Au at 1.5 A-GeV are rather well 
described in case of calculations without any medium effects as well as with all medium 
effects included. The dominant medium effect is a 'softening' of the antikaon spectra 
which is more in line with the data. Furthermore, the full in-medium calculations agree 
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very well with the rapidity spectra for semi-central collisions of Ni+Ni at 1.93 A-GeV 
while the spectra are underestimated by up to a factor of 2 for 'free' transition matrix 
elements and spectral functions. In addition, the centrality dependence of the K~ / K'^ 
ratio for Au + Au reactions at 1.5 A-GeV is rather well in line with the preliminary data 
of the KaoS collaboration for our full G-matrix calculations including pion dressing. The 
latter limit also provides the best description for the angular distributions in the cms, 
especially for non-central Au + Au reactions at 1.5 A-GeV (Fig. 18), which we interpret 
as a consequence of the lower amount of s-quark exchange reactions in the full G-matrix 
calculations. 

Nevertheless, in all limits considered in this work there is no convincing description 
of all spectra simultaneously. This failure might be cither attributed to larger systematic 
errors in the experimental data - as found for production in p -|- ^4 reactions from 
1-2.5 GeV [90] - or to an inadequacy of the many-body schemes adopted for the coupled- 
channel G- matrix calculations, especially the 'pion dressing'. Additionally, there is no 
direct experimental test of the hyperon dynamics incorporated in the off-shell transport 
approach since no explicit hyperon spectra are available for the systems measured at the 
SIS accelerator. These issues can only be settled by new experimental data in a wider 
kinematical regime with emphasis on low antikaon momenta in the center-of-mass frame 
since the in-medium effects are found to be most pronounced in this region as established 
experimentally also in production from p-\- A reactions for bombarding energies from 
1-2 GeV [91]. 
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Figure 1: The real part of the antikaon potential Re as a function of ^/s (Eq. (31)) for 
different nuclear densities and momenta = 0? 150, 300, and 500 MeV/c, respectively. 
The arrows show the pole mass of the antikaon in free space. 
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Figure 2: The imaginary part of the antikaon potential Im Ux as a function of ^/s 
(Eq. (31)) for different nuclear densities and momenta = 0, 150, 300, and 500 MeV/c, 
respectively. The arrows show the pole mass of the antikaon in free space. 
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Figure 3: The antikaon spectral function A{j)j^, y^) (Eq. (32)) for different nuclear den- 
sities and momenta -px = 0; 150, 300, and 500 MeV/c. 
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Figure 4: The K p cross sections to the final channels S vr , S^vr , S tt^ and Att 
according to the G-matrix calculations in 'free' space in comparison to the data from [62] . 
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Figure 5: The transition probability -Pi+2^3+4('S) (Eq. (34)) for the channel K~p K~p 
as a function of and the nuclear density p (in units of po) for a momentum Pk = 
relative to the nuclear matter rest frame. The l.h.s. corresponds to a calculation without 
pion dressing whereas pion dressing is included in the G-matrix calculations of the r.h.s. 
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Figure 6: The transition probability Pi+2^3+4(s) (Eq. (34)) for the channel K~p S°7r° 
as a function of ^/s and the nuclear density p for a momentum pj^ = relative to the 
nuclear matter rest frame. The l.h.s. corresponds to a calculation without pion dressing 
whereas pion dressing is included in the G-matrix calculations of the r.h.s. 



32 




Figure 7: The transition probability -Pi+2^3+4('5) (Eq. (34)) for the channel Avr'' —>■ K~p 
as a function of ^/s and the nuclear density p for a momentum pj^ = relative to the 
nuclear matter rest frame. The l.h.s. corresponds to a calculation without pion dressing 
whereas pion dressing is included in the G-matrix calculations of the r.h.s. 
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Figure 8: The transition probability Pi+2^3+4(s) (Eq. (34)) for the channel Avr^ Att^ 
as a function of a/s and the nuclear density p for a momentum pj^ = relative to the 
nuclear matter rest frame. The l.h.s. corresponds to a calculation without pion dressing 
whereas pion dressing is included in the G-matrix calculations of the r.h.s. 
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Figure 9: The experimental data for the reactions pp —>■ K^K^pp, pp K^K~pp (from 
Refs. [70, 62, 71]) and tt~p K~K^p (from Ref. [62]) in comparison to our parametriza- 
tions (sohd hues). 
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Figure 10: The in-medium cross sections for K' production from A^A^ and ir'p collisions 
for on-sliell nucleons and kaons in the final state, however, employing the antikaon spectral 
functions from the G-matrix approach in Section 3 as a function of ^/s for different nuclear 
densities ranging from 0.25 po to 2.25 po in comparison to the cross section in free space 
(solid lines). 
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Figure 11: The time evolution of the antikaon spectral distribution (Eq. (38)) in central 
Ni + Ni collisions at 1.8 A-GeV. The l.h.s. displays the antikaons stemming from the ttY 
channel whereas the r.h.s. shows the antikaons emerging from baryon-baryon collisions. 
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Figure 12: The differential inclusive (upper part) and K~ spectrum (lower part) for 
the system C + C at 1.8 A-GeV and 6cm = (90 ± 10)" in comparison to the data from Ref. 
[81]. The dashed line in the upper part reflects the result of the transport calculation 
after averaging the spectra over the angle Q in the center-of-mass system, while the 
solid line displays the calculated spectrum for 9cm = (90 ± 10)°). Lower part: The dashed 
line corresponds to a 'free' calculation, the solid line with open triangles to a (^-matrix 
calculation including pion dressing whereas the solid line with open circles results from a 
G-matrix calculation without pion dressing. 
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Figure 13: The differential inclusive (upper part) and K~ spectrum (lower part) for 
the system Ni + Ni at 1.8 A- GeV and 6iab = (44 ±4)° in comparison to the data from Refs. 
[82, 83]. The dashed line in the upper part reflects the result of the transport calculation 
after averaging the spectra over the angle Q in the center-of-mass system, while the 
sohd line displays the calculated spectrum for 9iab = (44 ± 4)°). Lower part: The dashed 
line corresponds to a 'free' calculation, the solid line with open triangles to a G-matrix 
calculation including pion dressing whereas the solid line with open circles results from a 
G-matrix calculation without pion dressing. 
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Figure 14: The differential inclusive K'^ (upper part) and K~ spectrum (lower part) for 
the system Au + Au at 1.5 A-GeV and 6cm = (90 ± 10)° in comparison to the preliminary 
data from Ref . [84] . The dashed line in the upper part reflects the result of the transport 
calculation after averaging the spectra over the angle Q in the center-of-mass system, 
while the solid line displays the calculated spectrum for 9cm = (90 ± 10)°). Lower part: 
The dashed line corresponds to a 'free' calculation, the solid line with open triangles to 
a G-matrix calculation including pion dressing whereas the solid line with open circles 
results from a G-matrix calculation without pion dressing. 



40 



10 



10 - 



10 



10 




5 -3 



■1.0 



-0.5 



0.0 



0.5 



1.0 



1 ' 1 ' 1 ' 1 ' 1 

: Ni+Ni->K +X, 1.93 A GeV, semi-central (b<4.5 fm) : 



J J'.... K 



r / rviicc 

rJl 

^ A fp^ — ° — ^ in-med. (without pi 

U f - ' 

±LJJ I 



pion dressing) 
K in-med. (with pion dressing) 

J 



■1.0 



-0.5 



0.0 

y 



0.5 




cm 



Figure 15: The rapidity spectrum of (upper part) and K~ mesons (lower part) for the 
system Ni + Ni at 1.93 A-GeV and semicentral colhsions (6 < 4.5 fm) in comparison to the 
data from Ref. [85]. The dashed hne (lower part) corresponds to a 'free' calculation, the 
solid line with open triangles to a G-matrix calculation including pion dressing whereas 
the solid line with open circles results from a G-matrix calculation without pion dressing. 
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Figure 16: The inclusive momentum spectrum of K~ mesons for the systems C + C 
Ni + Ni and Au + Au as a function of the antikaon momentum in the center-of-mass 
system pcm- The dashed hues correspond to a 'free' calculation, the solid lines with open 
triangles to G-matrix calculations including pion dressing whereas the solid lines with 
open circles result from G-matrix calculations without pion dressing. 
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Figure 17: The K~ /K^ ratio as a function of centrality, expressed in terms of the number 
of participating nucleons Apart, for the system Au + Au at 1.5 A-GeV and 0cm = (90 ±10)° 
including a cut in momentum pcm > 0.35 GeV/c. The prehminary experimental data have 
been taken from Ref. [84]. The dashed line corresponds to a 'free' calculation, the solid 
line with open triangles to a G-matrix calculation including pion dressing whereas the 
solid line with open circles results from a G-matrix calculation without pion dressing. 
The fluctuations of the lines are due to the limited statistics. 
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Figure 18: The (upper part) and angular distribution (lower part) for semi-central 
(l.h.s.) and non-central (r.li.s.) Au + Au collisions at 1.5 A-GeV. All angular distributions 
are normalized to unity for cos = 0. The preliminary experimental data have been 
taken from Ref. [88]. The dashed line corresponds to a 'free' calculation, the solid line 
with open triangles to a G-matrix calculation including pion dressing whereas the solid 
line with open circles results from a G-matrix calculation without pion dressing. The 
fluctuations of the lines for the K~ angular distributions are due to the hmited statistics 
again. 
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